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Abstract 

We present Molecular Dynamics simulations of the thermal glass transition in a 
dense model polymer liquid. We performed a comparative study of both constant 
volume and constant pressure cooling of the polymer melt. Great emphasis was laid 
on a careful equilibration of the dense polymer melt at all studied temperatures. 
Our model introduces competing length scales in the interaction to prevent any 
crystallisation tendency. In this first manuscript we analyse the structural properties 
as a function of temperature and the long time or a-relaxation behaviour as observed 
in the dynamic structure factor and the self-diffusion of the polymer chains. The 
a-relaxation can be consistently analysed in terms of the mode coupling theory 
(MCT) of the glass transition. The mode coupling critical temperature, T c , and the 
exponent, 7, defining the power law divergence of the a-relaxation timescale both 
depend on the thermodynamic ensemble employed in the simulation. 

1 Introduction 

Understanding the glass transition in supercooled materials is at the same time a great 
challenge in condensed matter theory [p], [2| and of high technological importance ||. 
Polymers constitute a class of materials with a very small crystallisation tendency. The 
ubiquity of amorphous polymeric materials has made them a longstanding focus of the 
experimental characterisation of the glass transition M as well as efforts to derive models 
of this transition ||. Traditionally this work focused on the temperature range in which 
the typical relaxation times in the material are macroscopic, i.e. in the range of seconds. 

With the development of the mode coupling theory of the glass transition 0, |6|, 0, [8|, |{J 
the interest was shifted to the temperature regime of the undercooled liquid and relaxation 
times in the ns to /xs range. This ignited a tremendous effort over the last decade IT, [12 



on the side of experiment and computer simulations to test the predictions and range of 



validity of this theory on all kinds of glass forming materials. Originally this theory has 
been developed for hard sphere liquids but it has been applied to and claimed to have 



been tested on as diverse materials as colloids |13|, p| , ionic glasses [[TjJ, molecular [TIJ 



liquids and polymers fT7 |. The emerging picture seems to be, that the theory can be 



applied in situations where the glass transition is determined by the repulsive part of 
the intermolecular interactions, i.e. where there are no site specific attractive interactions 
which could give rise for instance to network formation like in Si02- 

On the computer simulation side there has been a very detailed test of the MCT on 
a mixture of Lennard- Jones particles [19[ which belongs to the class of materials 
discussed above. The glass transition in polymer melts has been studied in great detail 



with Monte Carlo simulations of the bond fluctuation model |2T], |22|, [24], [25], [26], [27 
which can also be applied to the modelling of real polymers p8| . The combination of 
lattice model and Monte Carlo method necessarily means that one in general studies the 
glass transition at constant volume and that one has completely neglected inertia effects 
in the short time dynamics. Any motions on scales smaller than a lattice spacing are com- 
pletely eliminated, of course. Both of these drawbacks can be remedied by resorting to a 
Molecular Dynamics simulation of a continuum model. Work in this direction mainly used 



atomistic polyethylene-like models. Early work |29|, [30| focused on the glass transition as 



a phenomenon of macroscopic time scales, observing the break in the dependence of the 
specific volume on temperature. This study as well as later work pi , |32| used high quench- 



ing rates loosing one of the main advantages of such polymeric glass formers, namely the 
ability to equilibrate chain conformations and local packing in the amorphous state with- 
out intervening crystallisation tendency. In the series of works P0| , |5T| , no systematic 
study of quenching rate effects or the degree of equilibration at different temperatures 
was reported on. Especially for the dynamic structure factor it has been shown WA E3, 
however, that the observed behaviour depends strongly on the degree of equilibration one 
has achieved. 

Thus we decided to perform a systematic study of the glass transition in a polymer 
melt using a simple coarse-grained polymer model in the continuum consisting of Lennard- 
Jones particles connected by nonlinear springs ||33|| . This model can be equilibrated with 
respect to local packing and chain conformations down to temperatures well in the regime 
of the undercooled liquid. All our results on the dynamic properties of the polymer melt 
are therefore equilibrium dynamics. We also decided to simulate this model under constant 
volume (NVT) as well as constant pressure (NpT) conditions to study the difference and 
similarities of the glass transition in these ensembles. Since both cooling methods follow 
different paths in the state space of the model we expect quantitative differences between 
the observations at constant volume and constant pressure, albeit qualitatively similar 
behavior, as was seen experimentally p4}| . For a direct quantitative comparison one would 
need, for instance, a whole set of constant pressure cooling curves. 

The remainder of this paper is organised as follows. In Section |2| we will discuss our 
model and the simulation procedure for the NVT as well as NpT simulations. Section 
[3] will present a comparison of the static properties of the melts in the NVT and NpT 
simulations. In Section ^ we will discuss our results for the a-relaxation dynamics as 
observed in the dynamic structure factor and the self- diffusion of the polymer chains and 
Section [| will present our conclusions. 

A detailed analysis of the /3-relaxation regime predicted by mode coupling theory along 



the lines of Refs. |22|, |23| will be presented in a forthcoming work. 



2 Model and Simulation Technique 

For modelling the inter- and intramolecular forces we used a bead-spring model derived 



from the one suggested by Kremer and Grest |33[ and also used in several recent simu- 
lations 1 35 , 36]. We however included here also the attractive part of the Lennard- Jones 
potential, since previous work |24], |25] had shown that without such an attraction the 



model would produce a negative thermal expansion coefficient. 

Each chain consisted of 10 beads with mass m set to unity. Between all monomers there 
acted a truncated Lennard- Jones potential: 
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where C was a constant which guaranteed that the potential was continuous everywhere. 
Since it was not our aim to simulate a specific polymer we used Lennard- Jones units where 
e and a are set to unity. Note that this means that all quantities are dimensionless. In 
addition to the Lennard- Jones potential a FENE backbone potential was applied along 
the chain: 
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The parameters of the potential were set to k = 30 and Rq — 1.5 guaranteeing a certain 
stiffness of the bonds while avoiding high frequency modes (which would require a rather 
small time step for the integration) and chain crossing. Furthermore with these parameters 
we set the favoured bondlength to a value slightly lower than the length favoured by the 
Lennard- Jones potential. Thus we introduced two different incompatible length scales in 
our system, which should help to prevent the emergence of long range order at lower 
temperatures. 

All simulations in the NVT ensemble were performed using a Nose-Hoover thermostat 
P7| , to keep the temperature at the desired level. In this technique the model system 
is coupled to a heat bath which represents an additional degree of freedom represented 
by the variable £• The equations of motion are 
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where Fj is the total force acting on particle i due to the potentials described above and 
Q represents the mass of the heat bath, while g is the number of degrees of freedom. Note 
that ( fluctuates around zero and can thus become negative. The mass Q has to be chosen 



with great care [[39], |40| since otherwise one may not obtain a canonical distribution. If 
for example Q is very large, the kinetic energy and therefore the temperature starts to 
oscillate with an undesired large amplitude. Instead of a canonical distribution one then 
obtains a two-peaked distribution. In principle any problems could be avoided by using 
a chain of thermostats [fITJ but that would have worsened the computational effort and 
was thus discarded. For optimal results the intrinsic frequency of the heat bath should be 
approximately equal to the intrinsic frequency of the model system of which a theoretical 
estimate was obtained by calculating the frequency of a particle in a fee lattice subjected 
to Lennard- Jones potentials. The intrinsic frequency of the heat bath is given by p0| : 



Setting u s equal to the theoretically obtained frequency and rearranging this equation 
yields an expression for Q. Note that Q depends explicitly on the temperature and there- 
fore had to be adjusted for every simulation temperature. During all simulations no sus- 
picious behaviour due to the choice of Q was observed. We also performed several Monte 



Carlo simulations using both a continuum configurational bias method (CCB) |42|, |43 
and so called smart reptation which were carried out at the temperature T = 1.0 in order 
to check the validity of the Molecular Dynamics algorithm and to investigate whether 
this could be a potentially faster means for obtaining equilibrated configurations. For the 
dense melts we studied we found, however, that the fastest way to equilibrate the system 
was to use our standard MD algorithm. The measured static properties, as obtained in the 
MC simulations, were in good agreement with the measured static properties of the con- 
figurations produced with the Molecular Dynamics algorithm. Furthermore the obtained 
energy distributions were similar. To check what influence the Nose-Hoover thermostat 
has on the Newtonian dynamics we also carried out some simulations in the microcanoni- 
cal ensemble and compared the results to the results of the simulations with Nose- Hoover 
thermostat. Both methods lead to the same results, for example the velocity autocor- 
relation function of the two simulations were identical (Figure 1). This means that the 
thermostat only has a weak influence on the Newtonian dynamics although one is able to 
tune the temperature with it very effectively. 

Starting configurations were obtained using the method proposed by Kremer and Grest 



33], |36|. Before subjected to the Nose-Hoover thermostat each thus generated configura- 
tion was propagated in the microcanonical ensemble (Q = oo). At the beginning of this 
step the velocities were rescaled several times in order to come close to the desired tem- 
perature range. In the next equilibration step the thermostat was switched on and the 
system was propagated until the mean square displacement of the centers of mass of the 
polymer chains had reached several Rg, R g denoting the radius of gyration. At this time 
all measured correlators already had decayed to zero. 

In order to speed up computational efficiency we applied a linked cell scheme combined 
with a Verlet table ] 3]. Because of the Nose-Hoover thermostat it was not possible to 



use a Velocity Verlet algorithm, instead we used a Heun algorithm j45] with a time-step 
of dt = 0.002. 

All simulations in the NVT ensemble were performed using 95 polymer chains each 
consisting of 10 monomers. The volume was held constant at V = 1117. 65(p = 0.85). Since 
the density was the same for all temperatures it was not necessary to generate starting 



configurations at every temperature, as it was for the simulations in the NpT ensemble, 
but one could use the ones generated and equilibrated at another temperature and equi- 
librate them again at the new temperature. Simulations were performed at temperatures 
T = 0.35,0.38,0.4,0.45,0.5,0.6,0.7,1.0 and 2.0. For statistical reasons ten different con- 
figurations were simulated at each temperature. The equilibration of a configuration at 
the lowest temperature required 30 x 10 6 MD steps or almost two weeks of CPU time on 
a IBM Power PC for each configuration. 

Since in the simulations of the NpT ensemble we wanted to keep the average pressure 
at p — 1.0 at all temperatures, the situation differed from the one of the simulation of 
the NVT ensemble. In a first step we used a MD algorithm which also allowed for volume 
fluctuations of the system |39], [|(| to obtain the average density of the system at a cer- 
tain temperature. These runs lasted up to 5 x 10 6 MD-steps. Afterwards, in a procedure 
analogous to the one described above, we used the found density to generate starting con- 
figurations which we used for NVT simulations. Note that we performed the simulations 
themselves at constant volume, but this procedure made sure that the average pressure 
was constant (within five percent) at all temperatures. This was done because NVT sim- 
ulations are computationally more efficient, and also because we observed better stability 
for NVT than for NpT simulations. At almost all temperatures we simulated 120 poly- 
mer chains again each consisting of ten monomers. Furthermore at each temperature ten 
different configurations were simulated and simulations were performed at temperatures 
T = 0.48, 0.5, 0.52, 0.55, 0.6, 0.65, 0.7, 1.0 and 2.0. 



3 Static Properties 

In this section we will discuss the static properties of the melts as a function of tempera- 
ture. The thermodynamic paths for our cooling processes in the NVT and NpT ensemble 
are shown in Fig. 2. In the NVT ensemble we started our simulation at modestly high 
pressure at a high temperature. Upon cooling the pressure decreases and becomes neg- 
ative around T = 0.7. This negative pressure has consequences which we will discuss in 
detail when analysing the structure factor of the melt. In the NpT ensemble we keep the 
pressure at ambient value and adjust the density upon cooling. 



3.1 Chain Conformations 

Let us now first look at the chain conformations upon cooling. The Hamiltonian we chose 
has no intramolecular bond angle part and therefore there is no tendency of our chains 
to become stiffer at lower temperatures. Consequently the size of the chains varies very 
little in our whole temperature range (R g = 2.10 at T = 0.35, R g = 2.23 at T = 2.0 
in the NVT ensemble and R g = 2.09 at T = 0.48, R g = 2.23 at T = 2.0 in the NpT 
ensemble). In Fig. 3 we show the behaviour of the structure factor of the chains in the 
NVT ensemble for the lowest simulation temperature. Also included is a Debye function 
57 1 calculated with the independently measured radius of gyration. The good agreement 



with the simulation data shows that the chains remain Gaussian on the large scale over 
the whole temperature range. The stiffness parameter Cn was in the range of 1.51 to 1.56. 



3.2 Packing Behaviour 



The effect of the two competing length scales we introduced into our model can be nicely 
seen looking at the monomer-monomer pair correlation function shown in Fig. 4. First of 
all we want to note that the pair correlation function shows no long-range ordering even 
at the lowest temperature we studied. The nearest neighbour peak, which is just a diffuse 
peak around r = 1.0 at high temperatures, splits in two upon cooling. The first of the 
peaks is due to the preferred intramolecular distance or bond length b = 0.96. The second 
peak is the preferred nearest neighbour position in the minimum of the intermolecular 
Lennard- Jones interaction at r m j n = 2e. 

For the NVT ensemble this real space behaviour transforms into the structure factor of 
the melt shown in Fig. 5. The amorphous structure is here manifest in the amorphous halo 
around q = 6.9 , which contains both intramolecular and intermolecular nearest neighbour 
contributions. With decreasing temperature the short-range intermolecular order increases 
and since this is the larger of the two length scales contributing to the amorphous halo, 
its position shifts to smaller ^-values at first. At lower temperatures, however, this shift 
is reversed and the peak moves to higher g-values as would be expected for thermal 
contraction of the sample. In the same temperature range a small peak at very small 
q- values (q ~ 1.7) develops. Both effects result from a microvoid developing in the system 
because we are in a range of negative pressure where the system would like to contract 
into a dense melt expelling the free volume. The microvoid, which contains up to around 
five percent of the simulation volume, can be identified both by visual inspection and 
numerical analysis. The position of the small q peak is related to the typical diameter 
of a microvoid. For the NVT simulations we therefore have to keep in mind that at the 
lowest temperatures our system is no longer homogeneous but contains a small amount 
of internal surfaces. 

For the NpT simulations this effect is of course absent as can be seen in Fig 6. In this 
case the structure factor shows the behaviour also seen experimentally with the amorphous 
halo moving to larger g-values due to the increased density at lower temperatures. 



4 Dynamic Properties 

In this section we will look at the temperature dependence of the largest relaxation time 
in the melt. For simple glass forming liquids this is called the a-relaxation time. This is 
the timescale at which a particle breaks free of the cage of its nearest neighbours and large 
scale structural relaxation becomes possible. For polymers this then also is the timescale 
on which local conformational rearrangements start to occur. The largest relaxation time 
in polymers, however, is the time for the overall renewal of the chain conformation, which 
is a factor of iV 2 (iV being the number of monomers in a polymer chain) larger for chains 
following Rouse dynamics and a factor of N 3 for larger chains where reptation effects 



have to be taken into account |47| . The temperature dependence of this longest relaxation 
time is determined by the temperature dependence of the prefactor in these scaling laws, 
which is the timescale for local conformational changes, which, as discussed, is enslaved 
to the a-process of the structural relaxation. 



4.1 Structural Relaxation 



We will discuss the structural relaxation in terms of the incoherent intermediate dynamic 
structure factor 

^ ) = (^ e ^ (<K(?)) )' (7) 

where M stands for the total number of monomers in the melt. As can be seen in Fig. 
7, which shows F q (t) at the peak position of the static structure factor, the intermediate 
dynamic structure factor starts to exhibit a two step relaxation process when lowering the 
temperature, the so called (3- and a-processes. In this paper we focus on the behaviour of 
the long-time a-process and leave a detailed analysis of the /3-relaxation to a forthcoming 
publication. 

Comparing the behaviour of the NVT simulations (Fig. 7a) and the NpT simulations 
(Fig 7b) we see that the slowing down of the structural relaxation and the development 
of the two step process occur for higher temperatures in the NpT ensemble. 

The behaviour of F q (t) at the first minimum of the static structure factor is qual- 
itatively the same, with the plateau region occurring at a smaller value of F q (t). We 
empirically define an a- relaxation timescale by the requirement: 

F g (r a ) = 0.3 . (8) 

In the undercooled liquid close to the mode coupling critical temperature the time- 
temperature superposition principle is expected to hold. Figures 8a and 8b show that 
indeed we find a superposition of the a-relaxation behaviour for the NVT simulation in 
the region 0.35 < T < 0.45 and for the NpT simulations it the range 0.48 < T < 0.6. 

One generally analyzes the a-relaxation behaviour also by fitting the empirical KWW 
function to the data: 

f(t) = Ae-W (9) 

In the temperature range where we found the time-temperature superposition principle 
to hold we find: 

(3 = 0.56 ± 0.04(WT); /3 = 0.7 ± 0.08(NpT) (10) 

The error bars are mostly due to the effect that one can change (3 by almost 15 % by 
changing the time interval over which one tries to fit the data. When one tries to fit the 
behaviour at higher temperatures the resulting values for (3 increase approaching unity at 
high temperatures. 

For the a-relaxation timescale the mode coupling theory of the glass transition predicts 
a power law divergence: 

r a (x(T-T c )^ (11) 

Fig. 9 shows that we indeed observe this behaviour with T c and 7 depending on the 
thermodynamic ensemble. For the NVT simulations we obtain T c = 0.32 ± 0.01 and 
7 = 2.3 ± 0.2 and for the NpT simulations we have T c = 0.45 ± 0.01 and 7 = 1.95 ± 0.15. 
For our model we therefore equilibrated our system to temperatures within 10 % of T c . 
Note that near T c in the NpT ensemble the density is much higher than for our choice of 
density in the NVT ensemble and thus the large difference of T c is not unexpected. 



However it should be mentioned that it is in principle also possible to fit our data with 
the well known Vogel-Fulcher-Tammann (VFT) equation: 



r exp( 



E 
T — To 



(12) 



We obtain T = 0.215 ± 0.02 and E = 1.1 ± 0.1 for the NVT-simulations and T = 
0.34 ± 0.02 and E = 0.93 ± 0.1 for the NpT-simulations. Close to the mode coupling 
critical temperature the VFT curve is well approximated by a power law divergence with 
exactly the same critical temperature and exponent as obtained from an independent 
mode coupling fit. Very close to T c the VFT curve flattens in comparison with the ideal 
mode coupling fit. This is again in accord with what would be predicted by an extended 
mode coupling analysis |22|] taking into account structural decay via activated processes. 



It is also a behaviour typically seen in experiment [4S] and simulations [50 . 



4.2 Polymer Self- Diffusion 

The overall conformational relaxation of polymer molecules can be conveniently analyzed 
by looking at their self- diffusion behaviour . For this purpose one can look at: 



9i(t) = {(r 3 At)-r^(0)) 2 j, (13) 

which describes the mean square displacement of the inner monomers (j labels different 
polymer chains). The analogous quantity in the center of mass reference frame of chain j 
(r^ m (i) being the position of the center of mass of polymer j at time t) is: 

g 2 (t) = ((r|(t) - 4 m (t) - 4(0) + rUO)) 2 ) • (14) 
The mean square displacement of the center of mass itself is: 

9s{t) = ((rjjt) - rL(0)) 2 ) . (15) 

And finally the mean square displacement of monomers at the free ends of the chains and 
its analogous quantity in the center of mass reference frame are defined as: 

94t) = ((r^(t) - ri„ d (0)) 2 ) (16) 



96® = (ttnd(*) - ~ KnM + rUO)) 2 ) • (17) 

Fig. 11a shows g\ to #5 measured at T — 1.0 and Fig. lib the same quantities at T = 0.35. 
For gi(t) one can distinguish several regimes. For short times t < 0.1 one observes a 
ballistic regime which is followed by a subdiffusive regime and finally a free diffusion 
regime. Such a behaviour is typical for polymer systems and is predicted by many theories, 
e.g. the Rouse model. As can be seen in Fig. lib, the situation is a little bit different at 
lower temperatures. The ballistic regime is now followed by a plateau like regime which 
precedes the subdiffusive one. Such a plateau regime is typical for glass formers and a 
sign of the onset of the structural arrest of the system. The height of the plateau is 
closely related to the size of the cage a particle is trapped in. Another difference to high 



temperatures is that the subdiffusive regime stretches out far more in time and, therefore, 
the free diffusion limit is reached only after long simulation times. 
The subdiffusive regime can be fitted using: 



g l (t) = a 2 {W 1 ty\ (18) 

We obtain X\ = 0.62 ±0.02 for all simulated temperatures. g±(f) behaves similar to gi(t). 
Here the diffusive regime is preceded by a subdiffusive one as well, which can be fitted by: 

94 (t) = a 2 (W 4 t) x \ (19) 

Again the exponent is approximately the same at all temperatures. We find £4 = 0.67 ± 
0.03. 

If our model chains would exactly follow the Rouse predictions the local monomer 
mobilities W\ and W4 should be equal and x\ and x 4 should be equal to 0.5. It is however 
a general finding from simulations, that for chains as short as ours one generally observes 



a smeared out crossover behaviour [21] from short time to long time diffusion instead of 
the predicted Rouse exponent. At later times gi(t), g^{t) and g^it) all show the expected 
simple diffusive behaviour: 

9i (t) = 6Dt, (20) 

where the self diffusion constant is the same for all i. At lower temperatures, especially for 
gi(t) and (74 (t), it is very hard to distinguish this regime from the preceding subdiffusive 
one. 

In Figure 12 D(T), Wi(T) and W±(T) are plotted against the temperature. As one sees 
for lower temperatures all quantities follow a power law behaviour: 



D oc (T — T C ) J (21) 
Wi oc (T-T c )\ z = l,4 (22) 

The critical temperature and the exponent are, within range of error, the same as those 
for the a-relaxation timescale in the respective ensemble. This shows the coupling of the 
conformational relaxation and diffusion of the polymer chains to the local structural a- 
relaxation in the melt as discussed in the beginning of this section. Note that it is again 
possible to fit the data with the VFT-equation using the same To and E (within range of 
error) as obtained when fitting the a-relaxation times. 



5 Conclusions 

In this paper we have presented a Molecular Dynamics simulation of the thermal glass 
transition in dense polymer melts. We have studied this transition at constant density as 
well as constant pressure. Our model is a coarse-grained bead-spring model with nonlinear 
springs connecting monomers along a chain and Lennard- Jones interactions between all 
monomers. In order to introduce packing frustration into the model we chose incompatible 
length scales for intra- and intermolecular nearest neighbour distances. All our results were 
obtained on well equilibrated samples. 



We showed that the static structure factor of our chains of length iV = 10 could 
be well described by a Debye function at all temperatures. The size of the chains is 
mostly temperature independent as we introduced no temperature dependent stiffness 
into the model. The two incompatible length scales in the Hamiltonian can be seen in 
a split of the first neighbour peak of the monomer-monomer pair correlation function 
at low temperatures. For the constant volume simulation the density we chose led to 
negative pressure for T < 0.7. This instability led to the buildup of microvoids taking up 
approximately 5 % of the simulation volume at low temperatures. The observed negative 
pressure is an indication that the void formation process is not fully completed on the 
timescale of the simulation. 

In this work we analysed the glass transition in terms of the a-relaxation process. 
The divergence of the a-relaxation timescale could be very well described by the power 
law behaviour predicted by MCT. Critical temperatures differ substantially and power 
law exponents differ slightly between cooling at constant volume and cooling at constant 
pressure. The divergence of the a-relaxation timescale also leads to a divergence of the 
largest relaxation time in the system which is the Rouse time for the overall renewal of the 
chain conformations. This divergence can be observed looking at the mobility of monomers 
on intermediate length scales (as measured by the rate constants W\ and W4) and the 
center of mass self-diffusion coefficient of the chains. All these quantities follow power law 
singularities with values for the critical temperatures and exponents in nice agreement 
with the behaviour of the a-relaxation timescale. The divergence could be equally well 
described by a VFT fit to the data and it could be shown that close to T c the MCT power 
law singularity is a tangent approximation of the VFT curve. 

A detailed analysis of the /3-scaling regime predicted by mode coupling theory of the 
glass transition will be presented in a separate publication. 
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Figure 1: Velocity autocorrelation function with and without (dashed line) Nose- Hoover 
thermostat. The two lines are practically identical meaning that the thermostat has only 
a weak influence on the Newtonian dynamics. Note that because we were using Lennard- 
Jones units all quantities shown are dimensionless. 



1.10 




Figure 2: Thermodynamic paths for the cooling experiments we performed in the simula- 
tion. 




Figure 4: Pair correlation function for the NVT simulations and the range of indicated 
temperatures. Curves for lower temperatures are shifted upward. 
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Figure 5: Structure factor of the melt for a set of temperatures in the NVT simulations. 
Note the appearance of a small peak at low q-values around T=0.6 which is due to the 
emergence of a micro void in the system. 




Figure 6: Same as Figure 5 for the NpT simulations. 




Figure 7: Intermediate dynamic structure factor at the first maximum of the static struc- 
ture factor for the NVT simulations (a) and the same for the NpT simulations (b). 




Figure 8: Same as Figure 7 but with times scaled by the a-relaxation timescale for the 
range of temperatures where the time-temperature superposition holds, (a) NVT ensem- 
ble, (b) NpT ensemble. 




Figure 9: Critical behaviour of the a-relaxation timescale close to the respective critical 
temperatures, T c = 0.32 in the NVT ensemble and T c = 0.45 in the NpT ensemble, plotted 
against T — T c . Open diamonds are for the NVT ensemble and closed diamonds are for 
the NpT ensemble. 




Figure 10: r a as measured in the NVT ensemble. Also shown are best fits with the VFT- 
equation and the predictions of MCT. The temperature range is shifted by T c = 0.32 to 
show the similarity of the predictions of the two equations for our data close to T c . Note 
that very similar plots can be obtained for the quantities discussed in Fig. 12. 





Figure 11: Mean square displacements as a function of time for T = 1.0 (a) and T = 0.35 
(b). (?3 describes the center of mass of the chain, gi and are center and end monomers 
respectively and #2 and g§ are the corresponding displacements in the center of mass 
reference frame. 
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Figure 12: Critical behaviour of the self-diffusion coefficient and the rate constants for the 
NVT (a) and NpT simulations (b). 



